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The concept of mutually unbiased bases is studied for N pairs of continuous variables. To find mu- 
tually unbiased bases reduces, for specific states related to the Heisenberg-Weyl group, to a problem 
of symplectic geometry. Given a single pair of continuous variables, three mutually unbiased bases 
are identified while five such bases are exhibited for two pairs of continuous variables. For N = 2, 
the golden ratio occurs in the definition of these mutually unbiased bases suggesting the relevance of 
nimiber theory not only in the finite-dimensional setting. 
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Mutually unbiased (MU) bases of Hilbert spaces with 
finite dimension d (as defined by Eq. ((T) below) are a 
useful tool. If you want to experimentally determine 
the state of a quantum system, given only a limited sup- 
ply of identical copies, the optimal strategy is to perform 
measurements with respect to MU bases [1]. To pass a 
secret message to a second party, you could use quan- 
tum cryptography to establish a shared key, a procedure 
which relies on MU bases in the space [2, 3] or 
|[?]. Sending a physical system carrying a spin through 
a noisy environment, the effect of the interactions on the 
state of the spin might be modelled by a specific quan- 
tum channel, conveniently described in terms of MU 
bases [5]. Finally, if you happen to be captured by a 
mean king, you might be able to meet his challenge by 
knowing about entangled states and MU bases [6, 7]. 

Many of the ideas which underlie physical concepts 
defined for discrete variables, that is, in a Hilbert space 
of finite dimension, survive the transition from spin op- 
erators to position and momentum operators. Quantum 
key distribution and quantum teleportation [9], for 
example, possess coimterparts for continuous variables 
jlO] which act on an infinite-dimensional Hilbert space. 
It is thus natural to inquire into MU bases for continuous 
variables which, in fact, naturally occur in Feynman's 
path integral formulation of quantum mechanics [11]. 
The properties of MU bases in an infinite-dimensional 
space might also provide new insights into the existence 
of complete sets of MU bases in spaces of finite dimension 
not equal to the power of a prime. 

Let us recall the definition of MU bases in and 
some of their properties. Two orthonormal bases Bb — 
{|V''>b=i...d and Bb' - are called MU if 



K > 



b^b', 
b^b', 



(1) 



other basis. The value of the overlap k is not arbitrary 
but one derives from ^ that k = 1/ by using the com- 
pleteness of the basis Bb, say. 



Schwinger 11211 describes how to construct two MU 
bases from any orthonormal basis of C*. They are found 
to be the eigenbases of two operators U and V each 
shifting cyclically the elements of the other basis. These 
operators satisfy commutation relations of Heisenberg- 
Weyl type, UV = e^'^^/'^VU , describing finite transla- 
tions in a discrete phase space (Tsl] . This approach has 
been generalized in fl?], where it is shown that if one 
finds n unitaries each cyclically shifting the eigenbases 
of all other unitaries then these n bases are MU. 

The number of MU bases in is limited tod+1. Such 
complete sets of MU bases were constructed first in the 
case of d being a prime number 1 15] and subsequently 
for d being a power of a prime [1]. For composite di- 
mensions d = did2 ■ ■ ■ dk, the factors being (powers of) 
different primes, it is currently unknown whether com- 
plete sets of MU bases exist [16]. Interestingly, compos- 
ite dimensions are rare for small values of d but predom- 
inate for large d. While it is possible to construct three 
MU bases for any d > 2, numerical evidence for d=6 
(the smallest composite integer) suggests that no four 
MU bases exist [17] and that many of their subsets are 
missing as well [18]. 

Let us now turn to continuous variables p and q, with 
[q,p] = ih, acting on the Hilbert space £2(1^) of square- 
in tegrable functions on the real line. The (generalized) 
eigenstates of position and momentum \q),q G M, and 
\p),p E M, respectively, are known to satisfy 



m = 



1 



Jqp/h 



(2) 



since each state of one basis gives rise to the same prob- 
ability distribution when measured with respect to the 



Thus, a natural generalization of ^ for bases {|V^s)}s 



2 



of an infinite-dimensional Hilbert space takes the form 



S{s^s') if 6 = 5', 

fc>o a b^b', 



(3) 



where the (5-normalization of the states reflects the fact 
that the labels s, s' are continous. Consequently, the 
eigenstates of the position and momentum operators 
provide an example of MU bases with k = 1/ \/2Trh. The 
appearance of generalized eigenstates is inevitable, be- 
cause no normalizable state exists which has a non-zero 
overlap with all elements of a countable orthonormal 
basis. 

Is it possible to find three or more MU bases for one 
pair of continuous variables? The momentum basis Bp 
results from a rotation of the position basis Bq by an 
angle tt/2. Thus, a third MU basis might be given by 
^i} = {\Qi})}qt,eR' the eigenbasis of the operator = 
qcosd + psind with eigenvalue q^, d G (0, 27r). Us- 
ing Wigner functions, one finds that the modulus of the 
overlap between states of Bq and B^ is 



(4) 



Thus, no basis B^ with d e (0, n/2) combines with Bq 
and Bp to give a triple of MU bases. 

There is, however, a symmetric choice of operators 
which does provide three MU bases. Consider the bases 
B± = {\q±)}q±eVL where q± = gcos(27r/3) ± psin(27r/3), 
obtained from rotating the position basis by the angles 
±27r/3, respectively. One finds 



1 



27rri|sin(27r/3)| ' 



(5) 

so that the triple B+, B-, and Bq is MU with overlap 
k = 1/ VTrfiVS in Q. Comparing this result with (|2), we 
realize that, for continuous variables, the constant k in 
(|3) may take different values for different MU bases. 

In spite of l|4|, it is possible to complement Bq and Bp 
with a third basis resulting in an asymmetric triple of MU 
bases. Consider Bq-p consisting of the eigenstates of the 
operator q — p = \/2q^i^ which cannot be obtained from 
g by a rotation due to the factor V2- Nevertheless, one 
finds (as stated in iH) that 

I('zk-P)l' = l(9-Pb)l' = l(p|9)l' = ^, (6) 

providing us with an asymmetric triple of MU bases. 

We now develop a systematic approach to MU bases 
for N pairs of continuous variables residing in product 
states. For = 1, we will be able to explain the ob- 
servations above. For N > 2, we will derive geometric 
conditions which express whether product-state bases 
are MU or not. A set oifive MU bases will be found ex- 
plicitly for two continuous variables. Subsequently, we 



will formulate conditions to be MU for bases which do 
not have to consist of product states only. 
The Heisenberg-Weyl operator 



f(a) = exp[i(Pg-Qp)/ri] 



(7) 



which translates the position of a wavefunction by Q 
and boosts its momentum by P, will play a central role. 
We consider the generator of an infinitesimal transla- 
tion in the direction a* ~ {Q, P), using the notation: 



P(7 - = a* • j • X with j 



-1 

1 



(8) 



where x = [q^p)*" . Denote the eigenstates of by |a, a) 
where a identifies a particular family of states and a la- 
bels an element of this family. They satisfy 



a^a a, a) 



a e 



(9) 



forming complete and (5-orthonormal families of states 
Sa- Their position representations are given by 

, , , 1 



{q\a, a) 



^rP{q-a/Py/2hQ 



y/2TTh\Q\ 



(10) 



if both P and Q are non-zero [20]. The scalar product 
between states from bases with labels a and b (7^ a) is 
found to be 



|(b,/3|a,a)p = 



1 



27rn|b* - j - a] 



(11) 



It is crucial for the following that the right-hand-side 
of i fTTIl depends only on the modulus of the symplectic 
product of the vectors a and b, which is equal to the (un- 
signed) area of the parallelogram defined by these vec- 
tors. The particular class of states considered here thus 
picks up the symplectic structure related to the commu- 
tation relations [xa_, ib] = —iha* • j • b. Note that Eq. |(10)| 
is consistent with (Til l since one has \q) = \ep,a), with 
e^ = (0,l). 

We are now in a position to derive sufficient condi- 
tions to have MU bases for pairs of continuous vari- 
ables X„ = {qn,Pn), n ^ I .. .N, with [qn,Pn'] = ihSnn', 

each pair x„ acting on a copy of £2(IR)- 

In a first step, we restrict the candidates for MU bases 
to iV-fold tensor products of the states in l|9ll. 



N 



|a, a) 



|ai,ai; 



|aAr, aN) 



(12) 



which define a complete and (5-orthonormal basis Sg. 
Using l(TT|l , the modulus of the scalar product of |a, a) 
and \h,(3) is given by 



N 



(a,a|b,/3)p = ]J |(a„,a„|b„,/3„) 



N 



(2^fi)^ JjK.j.b„| 



(13) 
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which can be written as 

|(a,a|b,/3)p = (27r?i)-^|a*.jjvb|-i (14) 

where 

a = ai «)...®ajv , (15) 

etc. Thus, in order that some bases Bg, Sg, . . ., be MU, 
the unsigned symplectic products between any pairs of 
the vectors a, b, . . . must take one and the same value, 

|a*-jwb| = |b*-jAr.c| = ... = A->0, (16) 

reducing the search for MU bases of product form ((TTI l 
to the search of product vectors a, b, ... in satisfying 
((T6] |. Having found a solution {a, b, . . .} for some value 
of the constant K, one finds a solution for any other pos- 
itive K' by rescaling each vector with the factor ^ K/K'. 

What is the maximal number of vectors satisfying ((T6] l 
for N pairs of continuous variables? Lacking a general 
solution, we consider this problem of symplectic geom- 
etry in some detail for = 1 and N — 2. 

N = 1: The constraints ([16) now read |a* j b| = |b* j- 
c| = |c* • j • a| = fc > 0. In fact, only three vectors need to 
be written here since one can show that it is impossible to 
have a fourth vector d of symplectic product k with a, b 
and c satisfying these conditions. This does not exclude, 
however, the existence of four or more MU bases built 
from an entirely different set of states. 

Working out the unsigned symplectic product of the 
vectors (0,-1), (1, 0), and (1, 1) leads to fc = 1, correctly 
reproducing the asymmetric solution presented in 
Similarly, the set of unit vectors (0, —1) and (±V3/2, 1), 
which is invariant under three-fold rotations, describes 
the symmetric configuration lO, with fc = \/3/2. These 
apparently different solutions are, in fact, closely re- 
lated. Consider all real 2x2 matrices m with unit de- 
terminant which, under conjugation, leave the matrix j 
invariant up to a sign, 

m*jm = ±j; (17) 

we will call these matrices unsigned symplectic. They 
clearly form a group which consists of the union of all 
real symplectic 2x2 matrices, denoted by Sp(l, R), and 
all these matrices multiplied by the matrix j in (|8j which 
(is not symplectic but) satisfies iflTl l with the minus sign. 
Due to ((l7t , symplectic products a* • j • b remain in- 
variant up to a sign under transformations of the form 
a m • a. Using unsigned symplectic transformations, 
it becomes possible to map the triple of vectors (0, — 1), 
(1,0), and (1,1) into a configuration with three-fold ro- 
tational symmetry which is equivalent to the three MU 
bases in (|5]l, up to a non-unitary scaling transformation 
as described after Eq. (161 1. 

N — 2: MU bases correspond to sets of product vec- 
tors a = ai (g) a2, b = bi §5 b2, . . . , with equal unsigned 



symplectic products. We now exhibit ^'z^e vectors which 
satisfy (161 1 with K = 1, namely 

CHO' 

Here the number R is the golden ratio, i.e. the positive 
solution oi E? — i? + 1. Each coefficient of the five 
vectors is a sum of integer multiples of the numbers 
1 and R. Hence, we find that the coefficients are ele- 
ments of a number field given by a quadratic extension 
of the integers (just as the field of complex numbers is 
an extension of the real numbers where i, the solution 
of + 1 = 0, plays the same role as R). Thus the link 
between MU bases and number theory which pervades 
the finite-dimensional case (surveyed in, e.g. |21]) also 
exists for continous-variables. Interestingly, MU bases 
for multiple qubits |22] or qutrits [23, 24] must contain 
entangled states, contrary to what we find here. 

In a second step, we construct MU bases for N con- 
tinuous variables from states not limited to the tensor 
products (l2ll . To do so we introduce metaplectic opera- 
tors which represent linear canonical phase space trans- 
formations in Hilbert space. Explicitly, consider the 
transformation A' = M • A, with A = {qi, - ■ ■ ,pn) = 
(Qi p) € ^'^^ srid M being a symplectic matrix of size 
2N X 2N. Then there is a unitary operator Um such that 
the translation operators T(A)-each a product of N op- 
erators of the form (Til-transform according to 

?7mT(A) = f (M • A);7m , (19) 

defining the metaplectic Um- If symplectic transfor- 
mations are composed, M M' • M", then the cor- 
responding metaplectic operators are composed in the 
same manner: C/m = C^m'C^M"- 

The use of metaplectic operators has been implicit in 
our earlier discussion where we obtained a set of states 
I a, a), satisfying O, which are MU with respect to the 
position eigenstates | q) . We now show that these states 
can be obtained directly by application of a metaplectic 
operator. Expand (T9l l in A and consider the linear term 
to obtain Umxa = xm-aUm- First, let TV = 1 and choose 
the symplectic matrix m such that m • a = (0, 1)*, so 
Xm a ~ <?• The eigenfunctions of ia in (9]l are then gener- 
ated by I a, a) — Umlq)- The symplectic matrix satisfying 
m-(Q,P)* = (0,l)*is 

where /^j e M parametrises a shear along the line defin- 
ing the states |a, a). It affects the phase of {q\a, a), but 
not its magnitude. 
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In order to discuss a more general construction of MU 
bases (with iV > 1) we use a general expression |25] for a 
metaplectic operators which correspond to a symplectic 
matrix M of dimension 2N , 



exp(i9) 



^|det(M-I)|y (27r?i) 



dA 



exp 



—A* • N • A 

2h 



f(A); 
(21) 

here is a phase which need not concern us further, 
N = ij(M + I)(M-I)"^ is a symmetric matrix, J = j© 
. . . © j a block diagonal generalization of j in and the 
integration is over the 2N dimensions of phase space, 
dA = dQi dQ2 ■ ■ ■ dP^- The matrices M and N may be 
written using blocks of dimension N x N, 



Mqq Mgp 



■vr _ ( '^qg 

N N 

1 ^ pq -L ^ pp 



(22) 

Consider the action of Um on A^-fold products of posi- 
tion eigenstates, |q) = jgi) ® . . . ® |gAr). Using ||2T] | and 
l|22) l, we find that states C/mIq) ^ |M,q) are unbiased 
relative to the position eigenstates, i.e.. 



l(q'|M,q)p = 



1 



1 



(27rft)^ |det(M - I)dct(Npp)| 



(23) 



It follows from the composition property of metaplec- 
tic matrices that states |M, q) with different M are, and 
that the magnitude of their overlap can be calculated by 
composing the underlying symplectic matrices: 



|(M,q|M',q')P - \{q\U^'UM,\q')\' 

= |(q|(M-lM'),q')P: 



(24) 



where the final expression is evaluated using | |23] |. Thus, 
the problem of finding MU bases associated with meta- 
plectic operators can be solved by finding symplectic 
transformations such that the resulting expressions on 
the right-hand-side of (|23l l take the same values. This 
may allow for a much larger set of MU bases than l(T6] l. 

Our principal results are conditions for bases related 
by a metaplectic transformation to be MU, namely |(16] | 
(for which we found a solution (|T8] l) and more generally 
(|24j (as yet unexplored). To conclude we point out open 
questions. Even in the case of = 1, it is not know 
whether more than three MU bases exist. To have only 
three MU bases would be slightly surprising as the limit 
of d oo passing through prime dimensions suggests 
the existence of an unlimited number of MU bases. The 
result (|4j confirms this expectation in a restricted sense- 
any pair of bases and B^' is MU but with possibly dif- 
ferent values for the overlap. Future studies will reveal 
whether the pairwise unbiased bases , I? G (0,7r/2)are 
as useful as a complete set of MU bases. 

It is also unknown whether the bases Bq and Bp can 
be supplemented by a third MU basis qualitatively dif- 
ferent from the one presented in Let the state 



be a member of such a basis. The conditions KglV")! = 
KpIV")! — l/V27r/i imply that its expansion coefficients 
in the position and momentum basis are constant mul- 
tiples of phase factors exp[if{q)] and exp[ig{p)], respec- 
tively, related to each other by a Fourier transform. 



dq . 



(25) 



Thus, if the only pairs of functions {f{q),g{p)) solving 
this integral equation consist of quadratic polynomials, 
then there are no MU bases beyond the ones exhibited 
so far. Unfortunately, the entire set of its solutions is not 
known to us. 
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